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The melting and magnetic disordering of the skyrmion lattice in the quantum Hall system at fill- 
ing factor v ~ 1 are studied. A Berezinskii-Kosterlitz-Thouless renormalization group theory is 
employed to describe the coupled magnetic and translational degrees of freedom. The non-trivial 
magnetic properties of the skyrmion system stem from the in-plane components of the non-collinear 
magnetization in the vicinity of skyrmions, which are described by an antiferromagnetic XY model. 
In a Coulomb gas formulation the 'particles' are the topological defects of the XY model (vor- 
tices) and of the lattice (dislocations and disclinations) . The latter frustrate the antiferromagnetic 
order and acquire fractional vorticity in order to minimize their energy. We find a number of 
melting/disordering scenarios for various lattice types. While these results do not depend on a par- 
ticular model, we also consider a simple classical model for the skyrmion system. It results in a rich 
T — phase diagram. We propose that the triangular and square skyrmion lattices are generically 
separated by a centered rectangular phase in the quantum Hall system. 



PACS numbers: 73.40.Hm, 75.40. Cx, 75.70.Ak 
I. INTRODUCTION 

A. General Remarks and Motivation 

For nearly two decades the study of the quantum Hall 
effect has been one of th&,most productive fields of con- 
densed matter physicsJja Recently, quantum Hall sys- 
tems with additional degrees of freedom have received 
considerable attentionfllij In the simplest case this de- 
gree of freedom is the electron spin. Ideas developed for 
this system can be adapted for other multi-component 
quantum Hall systems such as coupled layers, wide quan- 
tum wells, and quantum wells in semiconductors with 
several degenerate conduction band minima.1121 

Here, we study the effect of the electron spins. We. 
are motivated bjz-.recent nuclear magnetic resonancet-3 
and specific heaiO measurements exhibiting interesting 
finite temperature spin physics. At the Landau level fill- 
ing factor v = l/m, where m is an odd integer, the 
ground state of the twfrt-dimensional electron gas is a 
strong ferromagnetJlSElO i.e., the electronic spins are 
completely aligned even in the limit of vanishing Zeeman 
coupling. Perhaps surprisingly, the effective Zeeman field 
in this system is rather small because of band structure 
effects. In the following we consider the case v ~ 1. The 
low-energy excitations of the system at v = 1 are spin 
waves gapped at the Zeeman energy. However, the quan- 
tum Hall ferromagnet also has topologically non-trivial 
excitations, which are (not strictly correctly) referred to 
as skyrmion£3c^3 analogy with the Skyrme model 
in nuclear physics. taEil They can be thermally created 
in pairs of vanishing topological charge, similar to vor- 
tices in two-dimensional superfluids. Skyrmions are in 
fact also present in conventional itinerant ferromagnets 
such as iron but do not seem to have any observable 



consequences at low temperatures. What makes them 
crucial for quantum Hall ferromagnets is that the vanish- 
ing diagonal conductivity a xx — together with the fi- 
nite Hall conductivity a xy = ve 2 /h makes thft-akyrmions 
carry a quantized electrical charge of ±i/e.ll3'ti3 As a re- 
sult, skyrmions (or antiskyrmions) are present even ii?4b*| 
ground state if we move slightly away from v = lE3€3 
(For our purposes skyrmions and antiskyrmions behave 
identically and we refer to both as "skyrmions.") 

If one dopes electrons or holes into the two-dimensional 
electron gas at v = 1, they enter the system as skyrmions 
with charge =Fe but with more than one flipped spin. 
This effect can be seen in measurements of the magne- 
tization as a function of filling factoro The result is a 
non-collinear ground state since the magnetization in the 
vicinity of the skyrmion centers has components perpen- 
dicular to the magnetic field, which have a vortex-like 
configuration. In a collinear magnet the SO (3) spin sym- 
metry is broken to a SO (2) symmetry with respect to 
rotations around the magnetic field direction. It has 
one Goldstone mode, which is gapped at the Zeeman 
energy in the presence of a magnetic field. In a non- 
collinear magnet the SO (2) symmetry is further broken 
and there are two Goldstone modes, only one of which 
is gapped. The other, gapless mode corresponds to ro- 
tation of the non-collinear spin configuration around the 
magnetic field axis. Thus non-collinearity leads to the 
appearance of a new low-energy SO (2) ~ U(l) degree 
of freedom. In the long-wavelength limit the orienta- 
tion of the in-plane components of the magnetization of 
a skyrmion can be described by a single U(l) phase factor 



or by the angle 



,i<t> 

Moving further away from v = 1, more and more 
skyrmions are present and their interaction becomes im- 
portant. The skyrmion interaction contains a repulsive, 
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long-range Coulomb part and a short-range contribution 
related to the U(l) degree of freedom. The latter term, 
which we here call the magnetic interaction, favors a, 
tiparallcl alignment of the U(l) "spins." Brey et al. 
recognized that the magnetic interaction could lead to a 
square lattice of skyrmions instead of the usual triangu- 
lar (hexagonal) lattice since the square lattice allows Neel 
order of the U(l) degree of freedom whereas magnetic 
order is frustrated on the triangular lattice. Since the 
magnetic interaction is of short range, the Coulomb in- 
teraction dominates for small skyrmion densities and one 
expects a triangular crystal. The U(l) degree of freedom 
is then frustrated, with neighboring angles <p differing by 
±120°. On the other hand, if the density is sufficiently 
high, the energy gained from Neel ordering on the square 
lattice may outweigh the lost Coulomb energy. Further 
lattice types may also be possible, e.g., a centered rect- 
angular lattice, i.e., a square lattice stretched along the 
(11) direction, with Neel order. In Sec. |l| we employ a 
simple classical model of the skyrmion system to illus- 
trate which lattice types and structural transitions may 
be expected. We find a surprisingly rich phase diagram 
for the classical ground state. The classical ground state 
has also been investigated by Rao et alE3 employing a 
non-linear sigma model and by Green et alBB, who also 
study the lattice dynamics, see Sec. || below. 

If the skyrmion positions were fixed to the ideal lattice 
at all temperatures, the long-wavelength physics, in par- 
ticular the critical properties, would be well described by 
an antiferromagnetic lattice XY model. We would then 
expectj-a. Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion,t3E3 which separates a low-temperature phase of 
bound pairs of logarithmically interacting vortices and 
antivortices from a high-temperature phase where large 
pairs are broken in the sense that their interaction is com- 
pletely screened. These broken pairs, which essentially 
consist of free vortices and antivortices, destroy quasi- 
long-range order. In the skyrmion lattice, however, the 
positions are not fixed and the lattice itself can melt. 

The critical properties of a two-dimensional lattice 
without any internal degree of freedom .have been success^ 
fully described by Nelson and HalperinE3 and by YoungE2l 
applying BKT theory to dislocations and disclinations 
of the lattice. Melting of the triangular and square 
lattices proceeds in two steps, both well described by 
BKT theory,c3 unless one or both of these transitions 
is preempted by a first-order melting transition. At the 
lower transition bound pairs of dislocations with oppo- 
site Burger's vector decouple, leading to a liquid crystal 
phase with short-range translational order but persisting 
quasi-long-range orientational order of nearest-neighbor 
bonds. Note that dislocations are, like vortices, topologi- 
cal defects with logarithmic bare interaction. For the tri- 
angular lattice, the liquid crystal phase is called hexatic 
because is shows quasi-long-range order with respect to a 
sixfold rotational symmetry, whereas for the-square lat- 
tice it is called tetratic (fourfold symmetry) .Ea However, 
the square lattice is unstable in a system with Coulomb 



repulsion as the only interaction. At the higher tran- 
sition temperature, pairs of disclinations, i.e., defects in 
the bond orientation field, unbind, leading to an isotropic 
fluid. The bare disclination interaction is confining but 
the presence of free dislocations above the lowec-jnclting 
temperature leads to a logarithmic interaction. Ej 

For the skyrmion lattice, the U(l) degree of freedom 
may not only stabilize the square lattice structure,^ 
we also expect the short-range magnetic interaction 
to be strongly affected by lattice deformations, lead- 
ing to the coupling of magnetic and lattice degrees 
of freedom. There are two aspects to this coupling: 
First, the low-energy collective modes (XY spin waves 
and lattice vibrations) may be coupled. From gen- 
eral argumentst^lHH the dispersion of the lattice vibra- 
tion mode, usually called magnetophonon mode, close to 
the Brillouin zone center is expected to have the form 
k 3 / 2 . Cote et alH have performed time-dependent 
Hartree-Fock calculations for the collective mode spec- 
trum of square skyrmion crystals. The authors indeed 
find two distinct low-energy branches. One is linear in 
k and is interpreted as the gapless XY spin wave mode, 
whereas the other has the fc 3 / 2 magnetophonon disper- 
sion. There is no sign of mixing of these two modes at 
small k. In Sec. [I C we briefly show that the classical 



skyrmion model reproduces these features. 

Second, despite the fact that the collective modes are 
largely decoupled, the topological excitations (vortices, 
dislocations, and disclinations) may be coupled, leading 
to an interplay of the magnetic BKT transition and the 
BKT melting transitions. That something non-trivial 
happens is easily seen from Fig. [|: A dislocation in a 
square lattice leads to a phase mismatch of ±7r in the 
U(l) degree of freedom since the nearest-neighbor cou- 
pling is antiferromagnetic. Naively one could expect that 
this mismatch along the dashed line in Fig. |l| leads to a 
linear, confining term in the interaction of a dislocation 
pair. However, the magnetization can relax so as to mini- 
mize the mismatch energy, as shown in Fig. 0(b). In this 
relaxed configuration, the dislocation has acquired half 
a U(l) vortex and the dislocation interaction is again 
logarithmic. The main objective of this paper is to illus- 
trate this point further and to explore its consequences 
for melting and magnetic disordering of various skyrmion 
lattice types. 

The vorticity acquired by dislocations can be derived 
usingpa. gauge theory of elasticity following ideas of Cardy 
et alx3, who have studied a soft square lattice of anti- 
ferromagnetically coupled Ising spins. Dislocations frus- 
trate t he a ntiferromagnetic order is this system as well. 
In Sec. IB we sketch the gauge theory for the case of a 
square skyrmion lattice. 

Experimentally, the situation is less clear. The pres- 
ence of skyrmions for v w 1 but v ^ 1 is firmly es- 
tablished by magnetization measurements by Barrett et 
aZo, using optically pumped nuclear magnetic resonance 
(NMR) techniques. Less certain, the rvery large nuclear 
relaxation rate T^ 1 seen in this regimeta is interpreted in 
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terms of the gapless XY magnon mode0 This mode cou- 
ples strongly to the nuclear spins because of its large S x,y 
components and its gaplessness. This coupling opens a 
channel for rapid spfu-lattice relaxation of nuclear spins. 
In these experimentsli-3 the skyrmions are probably usu- 
ally in a liquid state. Nevertheless the gapless XY mode 
is presumably-still present as an overdamped mode. 

Bayot et a/.lij find a strongly enhanced specific heat C 
for v w 1, which suggests strong coupling between elec- 
tronic and nuclear spins so that the large specific heat of 
the nuclear spins is in fact measured. Again, a plausi- 
ble coupling mechanism is provided by the gapless XY 
magnons. The filling .factor dependence of C is consis- 
tent with this picture. ta The temperature dependence of 
C shows a sharp peak at very low T. This peak may 
indicate a skyrmion lattice melting transition. Our so 
far quite speculative interpretation is the following: Ne- 
glecting the skyrmions for a moment, the Zeeman energy 
of nuclear spins within the quantum well containing the 
two-dimensional electron gas is (Knight) shifted because 
of their coupling to the polarized electron gas. Outside 
of the quantum well there is no such Knight shift and 
the mismatch in the Zeeman energy prevents the spins 
within and outside of the well from coming into thermal 
equilibrium. If skyrmions are present in a liquid state 
they move around, leading to motional narrowing and 
an averaged, but still finite, Knight shift within the well. 
In a lattice state, however, there are regions around the 
skyrmions where the electronic magnetization is perpen- 
dicular to the external field and the Knight shift van- 
ishes. The nuclear spins outside of the well can come 
into equilibrium with the nuclei in these regions. Hence, 
specific heat measurements suddenly see the nuclei out- 
side of the quantum well when a skyrmion lattice forms. 
Below the transition the apparent specific heat drops off 
again, which may indicate that the coupling is strong 
only in the vicinity of the transition where critical slow- 
ing down causes the electronic motion time scale to pass 
through thepNMR time scale. To our knowledge, these 
experimental^ are the only ones showing signs of a finite- 
temperature phase transition in a single-layer quantum 
Hall system. Recent experiments using resonant inelas- 
tic light scattering off double-layer systemsEJ also show 
signs of a finite-temperature transitionOEJ 

The objective of this paper is to illustrate several of the 
points raised above, in particular we wish to explore the 
consequences of the vorticity acquired by lattice defects 
for the melting and magnetic disordering transitions of 
skyrmion lattices. In Sec. HI we study these transitions 
for several possible lattice types. We introduce a gener- 
alized BKT renormalization group theory for a Coulomb 
gas with more than one species of particle carrying more 
than one kind of charge. The specific lattice types we 
study are motivated by the possible ground states of the 
simple classical skyrmion model of Sec. Q. It should be 
kept in mind, however, that the BKT theory does not 
depend on any particular model. 



B. Gauge Theory of Elasticity 
for the Skyrmion Lattice 

Following Cardy et aZ.El we here formulate a gauge 
theory for the magnetic and elastic energy of a square 
skyrmion lattice, i.e., a soft square lattice with an XY 
degree of freedom. Other skyrmion lattice types can be 



treated similarly. Let S(R) 



(R) 



be the XY spin 



field, where R is a lattice vector. We define a magnetic 
order parameter 



$(R) = S(R) e 



i7r(r x +r B )/a 



(1) 



where r(R) is the actual position of the skyrmion be- 
longing to the ideal lattice vector R and a is the lattice 
constant. This picture breaks down in the presence of free 
disclinations, i.e., above the upper melting transition, be- 
cause then the x and y components of the position vector 
r are no longer well-defined. The additional phase fac- 
tor turns the spins on one sublattice through tt, thereby 
mapping the antiferromagnet onto a ferromagnet. In the 
presence of dislocations this is not possible: The phase 
factor is no longer unique and ^(R) cannot be both con- 
tinuous and single-valued. Using a continuum notation, 
the exchange energy is 



(2) 



Next, we define two translational order parameters 

* a (R) = e a ^l a , (3) 
tf„(R) = e i27rr » /a . (4) 

Although r(R) is not continuous and single-valued in 
the presence of dislocations, the fields ^ x ,y are. Contin- 
uum elastic theory yields the elastic energy of a square 
lattice,E3 



(5) 



where summation over repeated indices is implied and 
1 / dr.; dr 



2 KdRi dR t 



(6) 



The last term in the elastic energy would be-absent for 
a triangular crystal or an isotropic medium.E3 Express- 
ing the Uij in terms of derivatives of the fields ^ x ,y and 
integrating by parts, the elastic energy becomes 

+ (m + A ) [(d x * x )*(d x * x ) + {d y %)*(d y %) 

+ + A + V) [y x (d x V x )*%(d y %) 

+ %(d y %yK(d x y x )]). (7) 
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The long-range Coulomb repulsion in the skyrmion lat- 
tice drives the Laiswi-coemcient A to infinity — the lattice 
is incompressible. Ej'E3 Thus only the first term on the 
right-hand side of Eq. M) is relevant (the other terms 
may introduce constraints on ^ x , y , which we ignore in 
the following since they do not affect our argument). 

A translation in the x (y) direction leads to phase fac- 
tors in ty x (\fy) and in A spin rotation leads to a 
phase factor in <E> alone. We can express the symmetries 
under translation and spin rotation as gauge symmetries: 
We introduce three two-component gauge fields A , A x , 
and A y and write the energy H = H cx + H e \ as 



H 



d 2 R 



a 2 [i 



iA Q - % -A x - ~A„ ) <i> 

a 2 n 



(8) 



^ |(V - iAjVj + l(V - iAy)* y 



This Hamiltonian is invariant under gauge transforma- 
tions with respect to any of the three gauge fields: 



A -> A + V0 O , * 



4>. 



(9) 



A x ^A x + V0 x , $^e^/ 2 $, V x ^e i6 *y x , (10) 
A y -» Ay + V6 y , $^ e ^/ 2 $, 9 V -► e ie " %. (11) 

These three transformations correspond to spin rotation, 
translation in the x direction, and translation in the y 
direction, respectively. The matter fields $, ^ Xl and ^ y 
in the Hamiltonian, Eq. (||), are only coupled through 
the gauge fields. 

We now discuss the topological defects in this theory. 
Magnetic vortices, i.e., vortices in are threaded by 
one flux quantum with regard to Ao. Dislocations cor- 
respond to topological defects in $ x or depending 
on the Burger's vector orientation. The elementary de- 
fect in, say, ^ x is a unit vortex. It is threaded by one 
flux quantum in Aa;. This seems to make the field <£> 
multi-valued since its phase changes by 7r if one moves 
around the vortex. However, if the XY spin part S in 
4> itself contains half a vortex (or antivortex) the order 
parameter <I> is again single-valued and continuous. This 
corresponds to ±1/2 flux quanta in Aq. Thus we reob- 



tain the result already discussed in Sec. [A: Dislocations 
each acquire ±1/2 magnetic vortex. 



II. CLASSICAL MODEL FOR 
THE SKYRMION LATTICE 

In the present section we formulate a simple classi- 
cal model for the interacting skyrmion system. We use 
this model to obtain (a) the classical ground state of 
the skyrmion lattice for a wide range of values of the 
skyrmion density and the magnetic interaction strength, 
and (b) the spectrum of low-lying collective excitations. 
This model should represent the physics of the real 



skyrmion lattice at least qualitatively, and even quanti- 
tatively at low density. This section is meant to illustrate 
some of the properties to be expected for the skyrmion 
lattice without introducing irrelevant technical compli- 
cations. Furthermore, we wish to motivate the choice of 



lattice types studied on a more general level in Sec. Ill 



A. Model 

The main idea is to take the correct limit of the two- 
skyrmion interaction at large distances and treat the 
skyrmion system as a classical gas of point particles hav- 
ing this interaction at all separations. We thus keep only 
the respective leading-order terms for large separations 
of both the interaction contribution independent of the 
XY degree of freedom and of the contribution depending 
on this degree of freedom, and we neglect three-, four-, 
etc. body interactions. This model should be valid at low 
skyrmion densities. 

We start from the .classical non-linear sigma model 
for the magnetization,t!3E3 which has &jf&\ successfully 
anplied to quantum Hall ferromagnetsliSlid Abolfath et 
alE3 have recently discussed the applicability of this clas- 
sical field theory and compared its predictions with mi- 
croscopic results. The magnetization is represented by 
a normalized three-component vector field m(r). The 
relevant terms in the Lagrangian read 



L = 



4tt£ 2 



J d 2 r A[m] • d t m - y J d 2 r (d i m ll )(d l m p ) 



P9*^bB f 2 



2e ./ r-r' 



(12) 



where A[m] is the vector potential of a—magnetic mo- 
nopole at the origin in spin (m) spacejij d t is a time 
derivative, and 



5p 



8- 



eijm ■ (dim x t9jin) 



(13) 



is the topological (Pontryagin) density. Greek indices 
always run over three values and latin ones over two. The 
first term is the usual Berry phase, and the other three 
stem from exchange, Zeeman, and Coulomb interaction, 
respectively. The Coulomb term reflects the fact that 
skyrmions carry electrical charge. 

In the absence of Zeeman and Coulomb interactions, 
the ground-slate solution for a single skyrmion is known 
analyticallyE3 It is scale invariant and for large distances 
r = |r| from the skyrmion center the in-plane components 
of m fall off as rrij oc rj/r 2 . Rotation around the z 
axis gives a different ground-state solution, reflecting the 
completely broken SO (3) symmetry. 

Switching on the Zeeman interaction, scale invariance 
is broken since the Zeeman term prefers a small skyrmion. 
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Far away from the skyrmion center we can expand the 
exchange and Zeeman terms in the Lagrangian, Eq. (12), 
up to second order iivthe small in-plane components nij 
of the magnetization,Lil 



E = d 



Y(dim j )(d l m j ) 



99* ^bB 



The resulting Euler-Lagrange equation 
pg*p B B 



didiirij 



has a vortex solution withEil 



r 3/2 



nij = 



(14) 

(15) 
(16) 



for j = 1,2. This expression is valid for r» 1/k, where 
k 2 = pg* p B B /2p s . Again the in-plane components can 
be rotated through any angle </>. By inserting m into 
Eq. (]l^) it is seen that the energy-density contributions 
of both the exchange and the Zeeman term fall off as 
r~ 1 e~ 2KT . Taking the Coulomb interaction into account, 
its leading-order contribution to the energy density also 
behaves like r _1 e _2Kr '. Thus all three energy contribu- 
tions are equally relevant at large r and Coulomb inter- 
action does not destroy the functional form of Eq. ( |l6| ) 
but does change the value of k. 

In calculating the interaction energy we assume that 
the two-skyrmion state with one skyrmion at the origin 
and the other at s is well described by 



mj(r) 



mf^{v)+R jk {ct>)ml 



single 



(r- s) 



(17) 



for j = 1, 2. The component ni3 is determined by |m| = 
1. This ansatz only gives errors of higher order in e~ KS 
for large separations. Here, 



COS( 

— sin ( 



sm < 
cos< 



(18) 



rotates the in-plane m components of one of the skyr- 
mions through an angle <f>. We find the interaction po- 
tential by inserting m into the potential energy part of 
Eq. ( |l2| ) and subtracting the energies of two isolated 
skyrmions. We are interested in the limiting form for 
large separations s. 

The exchange contribution to the interaction is the 
only one depending on the angle <fi. Using a multipole 
expansion and integrating over r we find, to leading or- 
der, the exchange contribution -E e xch oc cos 4> e~ KS /\/~s, 
where the coefficient of proportionality is positive. The 
contribution from the Zeeman term docs not depend on 
<f> and is exponentially small for large separations. We 
neglect it compared to the Coulomb interaction, below, 
since we only keep the leading 4> dependent and the lead- 
ing <p independent term. For v — 1 the leading contri- 
bution from the Coulomb interaction is e 2 /es, where e is 
the dielectric constant of the material. There are higher 
multipole terms, which fall off at least as 1/s 3 and are 
neglected compared to the 1/s term. 



B. Ground states 

As long as the skyrmion density is small, we expect the 
interaction to be dominated by the two-particle large- 
separation terms found above. Our approximation for 
the energy per skyrmion of a skyrmion lattice is 



E = E C 



E 



XY 



with 



Er 



e 

2~c 



E X Y = 



R^O 

gxY 



R 

cos(0 R - 



R#0 



R 



(19) 

(20) 
(21) 



where R runs over all skyrmion positions in the lat- 
tice except R = 0, Eq is an infinite constant from the 
Coulomb interaction with the neutralizing background, 
gxY > denotes the strength of the magnetic interac- 
tion, (f>n is the angle of rotation of the skyrmion at R, 
and ^xy = l//c is the range of the magnetic interaction. 

The long-range Coulomb interaction is not easy to sum 
over. The main idea of haw to make this summation well- 
behaved is due to EwalcO and has,.beeii successfully ap- 
plied to two-dimensional crystals:c3tl§E!J The lattice sum 
is split into a rapidly converging part and a long-range 
part, which is mapped onto a rapidly converging sum 
over the refiiotpcal lattice. Here, we quote a more gen- 
eral resultj£9t!J which will be useful later: If the R are 
summed over a two-dimensional Bravais lattice then 



-ik-(R+s) ^ 



R + s s ^ 



R 1 1 G 

+ v^J2 e ~ 4kR $ (H R + s| 2 ) 

R^O 

+ ^/n$(Trns 2 ) - -, 
s 



|G + k| 2 

Ann 



(22) 



where G are the reciprocal lattice vectors, n is the num- 
ber density, and Q(x) = ^/7rixerfc(y / a:) with the comple- 
mentary error function erfc.C 2 ! Equation ( p2[ ) only works 
for a Bravais lattice, lattices with a basis need special 
consideration. The magnetic structure is irrelevant here, 
since the Coulomb interaction does not depend on 

The simple sum over Coulomb interactions is obtained 
in the limit s — > 0, k — * 0, where the two sums on the 
right hand side can be cast into one, 



Er 



E 

R#0 



erfc(y / 7mi?) 2e 2 v / n 
R e ' 



(23) 



Now all lattice sums are rapidly converging and we can 
calculate the energy accurately. We write the energy per 
skyrmion in a dimensionless form, 
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F : E = V CTfc (^ Fr ) _ 2 



,-r//9 



with 



a ; 
/3 



e 2 n 1/4 /eg X Y, 



(24) 



(25) 
(26) 
(27) 



which arc all dimcnsionless. Note that (3 2 is the density in 
units derived from the range of the magnetic interaction 
and f3/a 2 = S,xY9xY e2 i s a measure of the relative 
strength of the magnetic interaction and does not depend 
on density. We expect this model to be quantitatively 
correct for small densities, (3 2 <C 1. 

The classical ground state for given a, (3 is determined 
by minimizing the energy (24). Thus we have to compare 
E for all reasonable two-dimensional lattice structures, 
taking the magnetic order into account. Besides the tri- 
angular lattice with frustrated antifetromagnetic order 
and the square lattice with Neel ordered we have also ob- 
tained ground state energies for the simple rectangular, 
the centered rectangular, and the oblique lattice^ thereby 
covering all two-dimensional Bravais lattices,c2l all with 
Neel order, and the honeycomb lattice, which is also bi- 
partite but is not a Bravais lattice. We cannot strictly 
exclude the possibility of more complicated ground states 
but have not found any other likely candidate. In the 
cases of the rectangular and oblique lattices, the lattice 
is characterized by one and two, respectively, continuous 
parameters in addition to its space group (Bravais lat- 
tice type). For example, the simple rectangular lattice 
has the anisotropy 77, defined as the ratio of the lattice 
constants in the (10) and (01) directions, as an additional 
parameter. To find the ground state, these parameters 
have to be optimized. 

Equation ( p3| ) is not applicable to the honeycomb lat- 
tice since it is not a Bravais lattice. However, its Coulomb 
energy Eq can be expresseiin terms of the triangular- 
lattice Coulomb energy EqE3 Taking the different den- 
sities into account, we find the dimcnsionless energy 
Eg = (1 + v / 3)/(2v / 2)-Bc- The Coulomb energies of 
the parameter-free lattices are Eq = —1.95013 for the 
square lattice, Eq = —1.96052 for the triangular lattice, 
and Eq = —1.89371 for the honeycomb lattice. The first 
two were also found in Ref. |4p| . 

We map out the ground-state phase diagram in Fig. |^ 
by following the various transition lines, i.e., lines of 
equal energy of two lattice types £3 We then discard lines 
that do not separate two regions with different ground 
states. The thin lines denote continuous transitions, 
whereas the heavy lines show first-order transitions .l3 
Recall that our approximation becomes doubtful for (3 2 ~ 
1, i.e., towards the right edge of the diagram. 



The phase diagram, Fig. g, is quite rich. For exam- 
ple, there is a region where the ground state is a honey- 
comb lattice. In its region of stability, it is even less frus- 
trated than the square lattice for our model interaction. 
In the upper left corner we find a very anisotropic ground 
state consisting of widely separated chains of skyrmions. 
Another interesting feature is the critical point on the 
square-simple rectangular line. Probably more relevant 
for real systems is the appearance of a centered rectangu- 
lar phase (a square lattice stretched along the diagonal) 
everywhere between the triangular and square lattices. 
It should be possible to experimentally see this two-step 
transition upon varrying v. Real skyrmion systems prob- 
ably live in the lower part of the phase diagram since the 
magnetic interaction cannot be made arbitrarily large in 
experiment. The parameter f3/a 2 can be increased by re- 
ducing the Zeeman interaction and thereby increasing the 
skyrmion size £xy- Experimentally, this can be done by 
applying hydrodynamic pressure. It is easier to increase 
the Zeeman interaction, reducing 13/ a 2 , by applying a 
in-plane magnetic field component. We roughly estimate 
that (3 jo? is smaller than unity in real systems. The 
transition lines show an upturn to larger magnetic inter- 
actions at the right edge of the phase diagram, Fig. ||. 
Although this may be an artifact of our approximation 
/? <C 1, it is interesting to note that a similar reentrance 
of a triangular phase is found in Ref. 



24 



The phase diagram is rather robust against changes 
in the exact form of the magnetic interaction. For ex- 
ample, the phase diagram for a simple exponential mag- 
netic interaction is qualitatively identical to Fig. [| This 
robustness indicates that the errors made by neglecting 
higher-order terms in the magnetic interaction are typi- 
cally small.,— - 

Rao et alB3 use a variational classical non-linear sigma 
model approach to find the classical ground states of 
the skyrmion lattice. They only consider the square 
and triangular lattices and consequently do not find the 
other phases, in particular the centered rectangular lat- 
tice. Their method has the advantage that the skyrmion 
size is optimized for given density and lattice type. The 
single skyrmion magnetization used in Ref. |24| does not 
approach the correct limit at large distances but accord- 
ing to the above argument this should not change the 
results qualitatively. However, at low density r-a, ferro- 
magnetically ordered triangular lattice is foundo which 
appears to be inconsistent with the large separation limit 
of the exchange interaction, Eq. (|2l|). 



C. Dynamics 

We now briefly turn to the low-energy collective ex- 
citations of the skyrmion lattice. As noted above, the 
U(l) degree of freedom leads to the appearance of a gap- 
less XY spin wave mode, whereas displacements of the 
skyrmions lead to a magnetophonon mode. The usual 
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ferromagnetic spin wave mode is gapped at the Zeeman 
energy. This mode is expected to mix with the XY mode, 
except at k = 0. This effect cannot be reproduced by the 
present model where the S z spin components are com- 
pletely integrated out. Thus the magnon dispersion is 
only reliable for the long- wavelength acoustical modes. 

We denote the displacements of skyrmions from their 
ground state positions by u = {111,112) and the devia- 
tion of the angle <fi from its ground state by uq. Then 
we expand the potential energy Eq. ( |l9| ) up to second 
order in u^. To describe dynamics we also have to know 
the leading time-derivative terms in the Lagrangian. The 
term for displacements-|pan be derived from the original 
Lagrangian, Eq. (O).cSEil In the limit of vanishing Lan- 
dau level mixing the skyrmion mass vanishes so that the 
Berry phase term is the only relevant one. One does not 
normally find a second-order time derivative term in spin 
dynamics, but Hartree-Fock calculationsEZl clearly show 
that Mo obtains a mass, or rather a moment of inertia /. 
This can be understood as arising from having integrated 
out all the short-wavelength spin fluctuations in order to 
obtain an effective action for the collective coordinate uq. 
There is also a Berry phase term associated with uq but 
it is a total time derivative and thus irrelevant at the 
classical level. The kinetic terms in the Lagrangian are 
thus 



T= T -J2 <(RK(R) + \ E ^<(R)^"(R) 



(28) 



R,n 



R,ri 



Here, R is a Bravais lattice vector of the magnetic lattice, 
the superscript n = 0, 1, . . . selects one skyrmion of the 
lattice basis, / is lh.ejnoment of inertia of it l (R), and 
the coefficient r\ iscJ'U r\ = e*B, where e* = ±e is the 
skyrmion charge. Deriving the Euler-Lagrange equations 
and making a plane- wave ansatz, 



we find the equations of motionEl 



(29) 







Ilu 2 A™ 



for \i = 
i-quie^A™ for fi ^ 



(30) 



with the dynamical matrix 

d 2 



-ik-R 



R 



E(R- 



c — c 



+ B,<f> n -<l> n +8 ) 



s=0,s =0 



(31) 



Here, E is the potential energy per skyrmion, Eq. (19), 
c™ denotes the position of skyrmion n within the unit 
cell, and 4> n is its XY angle. 

We can now address the question of mixing of magne- 
tophonon and XY magnon modes. Matrix elements mix- 
ing displacements and rotations of the XY angle (/i = 



and v ^ or vice versa) stem from the magnetic in- 
teraction alone and contain a first-order derivative of 
cos(</>™ — 0°), i.e., sin(0™ — <fi°), as a factor. Thus for any 
lattice with (j> n £ {0, ir} (square, rectangular, oblique, 
honeycomb) these matrix elements vanish. Of the lattices 
considered above, only the triangular shows any mixing 
of XY magnons and magnetophonons in our model. 

Here we only show results for the square lattice. 
The contribution from the magnetic interaction can be 
summed directly, whereas for the Coulomb interaction 
we need the second derivatives of Eq. (|2^). The magne- 
tophonon and XY magnon dispersions for /3 2 = 0.1 and 
j3 jo? — 1 are shown along directions of high symmetry in 
the magnetic Brillouin zone in Fig. |^. There are only two 
instead of four magnetophonon modes in the magnetic 
(folded-back) zone since it™ and u% are canonically con- 
jugate and thus do not lead to independent modes. The 
dispersion of the lower magnetophonon mode for small 
k is indeed of the form u oc fc 3 / 2 , whereas the magnon 
mode is linear. Note that the optical magnetophonon 
branch shows a minimum at the zone centec-pa, feature 
previously seen in Hartree-Fock calculations.E3c3 We see 
that our simple model reproduces the main qualitative 
features of the collective mode spectrum and conclude 
that it captures the essential physics. 

The dispersion relations can be used to reintroduce at 
least part of the quantum effects into our model by taking 
the zero-point energy hw/2 for all modes into account .c2l 
This additional energy may favor certain lattice types 
and thus shift the transition lines. It may also lead to 
quantum melting. Furthermore, we could extract low- 
temperature thermodynamic properties, such as the free 
energy, and with its help study structural phase transi- 
tions at finite temperatures. The harmonic approxima- 
tion is, in principle, inconsistent with the description of 
melting, be it quantum or thermal, which is intrinsically 
non-linear. Nevertheless, estimates of melting temper- 
atures could be found by calculating the amplitude of 
vibrations in this approximation and using a Lindemann- 
type criterion.cil In particular, we expect quantum melt- 
ing due to soft modes in the vicinity of the continuous 
transition lines, e.g., the one between square and cen- 
tered rectangular lattices. 

The magnon and magnetophnnon dispersions have also 
been studied by Green et al.tB This work is not eas- 
ily comparable to real systems since the authors assume 
that the skyrmion shape is little affected by Zeeman and 
Coulomb interactions, which yields an incorrect expres- 
sion for the interaction at large distances. Furthermore, 
they introduce a mass term into the displacement equa- 
tion of motion and their magnetophonon frequencies are 
inversely proportional to the mass, which seems doubt- 
ful since the mass vanishes for vanishing Landau level 
mixing so that the physics is dominated by the Berry 
phase term. Two lattice types are considered in Ref. |2^: 
what we call the square lattice and the centered rectan- 
gular lattice, the latter with the angle between primitive 
lattice vectors fixed to 7r/3. The latter is structurally 
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identical to the triangular lattice but has only two mag- 
netic sublattices. The authpss dispute the existence of a 
frustrated triangular phased but since they do not con- 
sider a triangular lattice with a three-skyrmion basis it 
is clear that they cannot find it. 

If we go to higher temperatures, the harmonic approx- 
imation breaks down. In particular, topological excita- 
tions become important. They are believed to lead to the 
ultimate melting of the skyrmion crystal. In the following 
section we describe possible scenarios of melting. 



III. BEREZINSKII-KOSTERLITZ-THOULESS 
THEORY FOR SKYRMION LATTICE MELTING 

In the present section, the central one of this pa- 
per, we discuss the melting transitions and intermediate 
phases of the lattice types discussed above. We use the 
framework of a suitably, generalized BKT renormaliza- 
tion group theoryo'LJo since this theory is known to 
work well for the simpler problems of (i) an XY model 
on a rigid lattice and (.ii) a soft lattice without additional 
degrees of freedom.ESEj More specifically, we employ a 
Coulomb gas language. BKT theory is a static theory 
so that the unconventional kinetic terms in the skyrmion 
Lagrangian, see, e.g., Eq. (p8|), do not affect the results. 
We stress that this approach does not make reference to 
any special model for the interacting skyrmion system. 
There is one important caveat: The theory of Refs. |28| 
and ^9] describes the system well if it shows BKT melt- 
ing but does not say whether it actually does. The upper 
or both BKT melting transitions may be preempted by 
a first-order transition. 



A. Multiple-charge Coulomb gas 

We first introduce a general model, which contains all 
relevant skyrmion lattices as special cases. This model 
is a two-dimensional continuum Coulomb gas with more 
than one species of particles carrying more than one 
charge. The dislocations and vortices are treated as clas- 
sical point particles with logarithmic interactions. From 
continuum elasticity theory on finds that for the triangu- 
lar lattice the interaction between two dislocations with 
Burger's vectors, h^-^iad b2 and separation vector r is 
proportional tceZH 



|r| (b^-rXba-r) 

bi • b 2 In 1 - . 

r v z 



(32) 



where the length scale r is given by the lattice spac- 
ing. For less symmetric lattice types this expression is 
not strictly correct but the leading, logarithmic term 
is always present rB We here only keep the logarithmic 
term. The qualitative behavior and the universal u 
in the stiffness are known to be unaffected by this! 28 



[However, the temperature dependence of the correla- 
tion length above the dislocation unbinding transition 
changes for the triangular lattice as a result of both 
the sub-leading term in Eq. ( |32"|) and the appearance 
of triple±a_pf dislocations with vanishing total Burger's 
vector.BEfl 

As discussed in the introduction, dislocations attract 
partial vortices to minimize the energy resulting from the 
mismatch in the antiferromagnetic order. In a Coulomb 
gas language we have three charges: the vortex strength 
and the x and y components of the dislocation Burger's 
vector. These three charges correspond to the three 
gauge fields discussed in Sec. IB. A similar description 



can be used for the possible upper melting transition from 
a liquid crystal to an isotropic fluid. This transition is 
thought to be due to unbinding pairs of disclinations, 
which may again bind partial vortices. 

The model is defined as follows: There are N species 
of particles, counted by 11 = 1,...,N, which carry 
M charges g^,...,^. Each particle has an antipar- 
ticle with all charges inverted, = —q™, where we 
use the notation n for the species of the antiparticle. 
The charges interact via the two-dimensional logarithmic 
Coulomb potential. The interaction between two parti- 
cles of species n and n' at positions r and r' is then 



M 



V 



(33) 



The charges have units of a square root of energy — the 
strength of the interaction is contained in q™. Further- 
more, we assume (i) that only a particle and its antipar- 
ticle can form a pair that is neutral with respect to all 
charges and (ii) that the whole system is neutral with 
respect to all charges. The number of n particles and n 
antiparticles is then equal. Restriction (i) is in fact not 
crucial but simplifies the argumentation. 

The usual Coulomb gas model for vortices in the XY 
model or a superffuid film is the case N = 1, M = 1. 
Melting of a square lattice without additional degrees of 
freedom can be described by an N = 2, M = 2 Coulomb 
gas, which reduces to two independent, identical N = 1, 
M = 1 models since dislocations with Burger's vectors 
along the x and, y axes, respectively, do not interact. 
Tupitsyn et alE3 have considered an N = 2, M = 2 
model for merons in a double-layer quantum Hall system, 
but with a 1/r interaction for one of the charges. 

For the square skyrmion lattice we have M = 3 charges 
and N = 5 particle species. The charges correspond to 
the vortex strength (<?„), the x component of the Burger's 
vector (q 2 ), and its y component (q^)- The charges of the 
particles are given in table Gj. This table reflects the fact 
that dislocations bind +1/2 or —1/2 vortex. 

We, apply Kosterlitz' renormalization group formula- 
tionB of BKT theory to the general N, M model. This 
approach has the advantage of being more rigorous and 
mathematically more transparent than, e.g., the original 
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self-consistent screening approach.c3c3 The renormaliza- 
tion procedure of Ref. [27] consists of two steps: First, 
the smallest neutral pair is integrated out, and then the 
length scale r, the size of the smallest pairs, is rescaled. 
The procedure for the multiple-charge Coulomb gas is 
similar to the original caseEH and is not given in detail. 

The main approximation of the original theoryEH is 
that the two particles with smallest separation are as- 
sumed to always form a neutral pair. This is reason- 
able since two particles which do not form a neutral pair 
have the same charge and thus repel each other. In our 
case this is not generally true: two particles which are 
not a neutral pair can even attract each other and form 
a non-neutral bound state. However, this arrangement 
will attract other particles until it is totally neutral. The 
approximation that the smallest pair is neutral is thus 
equivalent to neglecting neutral arrangements of more 
than two particles. The same problem arises for the tri- 
angular lattice without XY degree of freedom, where 
three elementary dislocations can have vanishing total 
Burger's vector O Depending on the g™ there can be neu- 
tral triplets, e.g., (123) in the above example, whereas in 
the N = 1 case the simplest neutral arrangements except 
for pairs are quartets. Therefore the first neglected term 
in the renormalization group equations is of fourth or- 
der in the patticle fugacities for both the XY model on 
a rigid latticellj and the melting square lattice without 
XY degree of freedom, whereas it is of third order for 
the square skyrmion lattice. 

The grand canonical partition function of the multiple- 
charge Coulomb gas is 



JV 

*= e n 

A/i,...,A/jv »»=1 



(MJ) S 



27V„ 



d 2 n 



Dt(r) 



[ d\ 2M exp I +| £ f) C< In 

(34) 

where 7V n is the number of particles of species n (there 
is the same number of n antiparticles), TV = 53 n 7V n , the 
y n are fugacities, the ranges of integration Di(r) com- 
prise the whole plane but exclude configurations with 
two particles closer than r, and the double sum 
runs over all 27V particles and antiparticles. Pairs of size 
in [r, t + dr) are integrated out according to 



n / 



n / 

i=l J Di(T+dr) 



d 2 r k 



d 2 r„ 



i7 t jk7H<j JD h (r+dr) JD> 
I (I 1 !-, 8„ 77 . 



(35) 

/r<|ri — rj |<T-t-G0T 

Here, D' consists of the whole plane except for disks of 
radius r centered at all particles k ^ The only ap- 



proximation here is contained in the symbol S ni ^ . , which 
states that only neutral pairs are integrated out. Apply- 
ing this prescription to Eq. ( |34"| ) and rescaling r we obtain 
(cf. Ref. ||) 



z^ Zo £ n 

Mi,—,Nn »=1 



(7V„!) 2 \{r + d T f 



2-/V„ 



M 



l + [2-/3/2^(C) 2 ]^ 

v. m=\ I 



2M n -, 



2Af 



d Ti exp 



i=1 JD,(r+<lr) 



E 



71/ 



E c< 



m— 1 



2n 2 j2vi^ ( E ccYe c<) 

71=1 \rn—l / \m— 1 / 



1 \ r i- r ]\ 

x In ■ 



dr 



(36) 



Except for the irrelevant constant Zq, this partition func- 
tion is identical to the original one if we replace 



m — 1 



(37) 



M 



M 



N 



E « - 2tt 2 E yi-P 

m— 1 n" — 1 

x (E«')(E^'^) -■ ( 38 ) 

\m=l / Vm=l / 



If we define a (symmetric) stiffness tensor 



2tt 



(39) 



and express the scaling relations for y n and K nn i by dif- 
ferential equations, we obtain the generalized renormal- 
ization group equations 



dK„ 



N 



df 



-An 3 E yn»Knn»Kn"n', 



(40) 
(41) 



where I = lnr/r is the logarithmic length scale. The 
initial conditions for these equations are 



yl{i = Q) = Cle~ 2 ^\ 

/3£ ro C (o)^(o) 



2n 



(42) 
(43) 



where C n are constants of the order of unity and £^ orc are 
the core energies of one n particle. Note that the general- 
ized Kosterlitz equations and the initial conditions reduce 
to the standard BKT expressions for N = 1, M = 1. 







B. Skyrmion lattices 

We first consider the square skyrmion lattice (N = 5, 
M = 3) . At first glance this problem looks rather compli- 
cated, since it involves coupled differential equations in 5 
fugacities and 15 stiffness constants, taking the symme- 
try of K nn i into account. However, we can simplify the 
problem considerably by looking for further symmetries. 
In a first step we consider the symmetries of y n and K nn i 
at the minimum length scale £ = and check which of 
these survive for the renormalized quantities at £ > 0. 
From table | and the observation that the core energies 
of all species of dislocations should be equal, we see that 
at £ = there are at most seven independent quantities, 



Inserting this ansatz into Eqs. 
pendent equations, 



t) we obtain four inde- 



Vv = 


2/1, 






yd = 


2/2 = 


2/3 = 2/4 


= 2/5, 


K v = 


K n 






K d = 


K 22 


= ^33 = 


if 44 = if 55, 


K vd = 


if 12 


= -K 13 


= if 14 - -if 15 


K dd = 


K 23 


= ^45, 




Kd = 


K 2i 


= -K 25 


= -if 34 - if 35 



(44) 



The equations for these quantities can be read off from 
Eqs. d40| ) and (|4l|). It is easy to see that two quantities 
that are equal by symmetry at £ = remain equal to 
each other as we integrate away from £ = 0. Hence, the 
problem can be reduced to seven coupled equations, 

dy 2 Jd£ = 2(2- nK v )yl 
dyj/di = 2(2 -nK d )yl 
dK v /d£ = -^ 3 y 2 v K 2 v - 16^if 2 d , 
dK d /d£ = -4n 3 y 2 v K 2 d - 4n 3 y 2 d (K 2 d + K 2 dd + 2K' 2 d \ 
dK vd /d£ = -^ 3 y 2 v K v K vd (45) 

- ^ 3 y 2 (K d K vd - K dd K vd + 2K' dd K vd ), 
dK dd /d£ = +Air 3 y 2 v K 2 d - ^ 3 y 2 (K d K dd - if^ 2 d ), 
dK'Jdl = ~4n 3 y 2 v K 2 d - 8n 3 y 2 (K d K' dd - K dd K' dd ). 

In a second step we make an ansatz for the remaining if 
to reduce the number of independent quantities further. 
We guess that the renormalization of the interactions can 
be described in terms of independent screening of the two 
charges q v and q d alone (we will see below that this as- 
sumption is not correct for all lattice types). If it were 
true there would be only two independent stiffness con- 
stants J v and J d , which we choose so that 



■7„(0) 



Pi 2 M JM _ 0«2(o) 



(46) 



2tt ' " ' 2tt 

Since our ansatz has to work at £ = we have 

K = J 

K d = J d + J v /4, 

K vd = J v /2, (47) 
K dd = J d — J„/4, 
K' dd = J v /4. 



dy 2 Jd£ = 2(2- it J v )y 2 v , 
dy 2 d /dt = 2(2 -wJ d -nJ v /4)y 2 d , 
dJ v /d£ = -4TT 3 y 2 v J 2 - 4n 3 y 2 J 2 , 
dJ d /d£ = -87r 3 y 2 J 2 , 



(48) 
(49) 
(50) 
(51) 



and three that are linear combinations of these. Thus 
our ansatz is indeed correct. Of course, we could have 
made this ansatz directly for the if„„' without the step 
in between. We repeat that the leading neglected terms 
in Eqs. (48)- (|5l|) are of third order in y v . d . 

From Eqs. (|5C|) and (|5l]) we see immediately that if 
the dislocations proliferate, lim^oo y d = 00, both stiff- 
ness constants J v and J d go to zero for £ — > 00. This 
result reflects the fact that free dislocations do not only 
screen the dislocation interaction but also the vortex in- 
teraction since they carry vorticity. On the other hand, 
free vortices (lim^oo y 2 = 00) only lead to J v — > since 
vortices do not have a non-zero Burger's vector, and, con- 
sequently, Eq. ( |5l| ) does not contain y 2 . We now discuss 
the possible scenarios. 

Decoupled transitions. Fig. |4|(a): We start from low 
temperatures. At some temperature T v , vortices unbind 
and J v (£ —> 00) shows a universal jump from 2/tt to zero. 
(We now omit the £ argument when we refer to the limit 
£ — > 00.) At the same temperature the effective stiffness 
of the dislocation interaction, J d + J v /4, also shows a 
jump but at the high-temperature side of the jump 
we still have J d + J v /A — J d > 2/tt. There is no jump in 
J d alone. Above T v the system has only short-range mag- 
netic order but still quasi- long-range translational order. 
At some higher temperature T d the dislocations unbind 
and J d jumps from 2/7p-to zero. Both transitions show 
BKT finite size scaling,El J Vtd (£, T V}d ) 2/tt[1 + 1/(21)]. 

Vortex driven simultaneous transitions, Fig. ^(b): At 
T v vortices unbind: J v jumps from 2/n to zero. At T~ , 
the effective dislocation stiffness is J d +J v /4 > 2/tt, while 
J d alone is smaller than 2/n. With J v jumping to zero, 
we appear to have J d + J v /4 = J d < 2/tt at X!+, sec 
the cross in Fig. f|(b). However, from Eq. ( p9| ) we see 
that then lim^oo y d — 00, dislocations also proliferate, 
and J d jumps to zero. The physical reason is that with 
the vortex interaction screened the remaining dislocation 
interaction is suddenly too weak to bind dislocation pairs 
so that T d = T v . We espect BKT scaling only in J v . 

Dislocation driven simultaneous transitions, Fig. ^(c): 
At T d dislocations unbind: the effective stiffness J d +J v /4: 
jumps from 2/tt to zero. The vortex stiffness is still large 
at T7 , J v (T d ) > 2/n so that vortices would unbind at 
a higher temperature. However, since J d + J v /4 van- 
ishes at T d , so does J v so that T v = T d . Physically, the 
vortex interaction is suddenly screened since the prolif- 
erating dislocations carry vorticity. For this reason the 
magnetic transition can never take place at a higher tem- 
perature than the melting. There can be lattice order 
without magnetic order but not vice versa. We expect 
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BKT scaling in Jd + J v /4 and also find similar scaling in 
J v and Jd separately but with a non-universal value of 
J v>d (T d ) ? 2/tt. 

Above the dislocation unbinding transition at Td the 
system still shows orientational quasi-long-range order, 
whereas translational and magnetic order are of short 
range. This— .is the tetratic phase mentioned in the 
introduction.^!! It is characterized by free dislocations, 
which are in fact bound pairs of disclinationsHa For 
the square lattice elementary bare disclinations do not 
carry vorticity, as can be seen from Fig. ||. Disclina- 
tions dress with-free dislocations, leading to a logarith- 
mic interactionra This screening cloud of dislocations 
is expected to have vanishing total Burger's vector and 
vanishing total vorticity in order to minimize its energy. 
Thus dressed disclinations still have zero Burger's vec- 
tor and vortex strength and the disclination unbinding 
transition corresponds to an N = 1, M = 1 model with 
the dressed fivefold disclination as the particle and the 
dressed threefold disclination as its antiparticle. At a 
temperature Tdi SC > T d disclination pairs unbind and the 
system becomes an isotropic fluid. 

We now turn to the other lattice types. The triangular 
lattice has dislocations with Burger's vectors along any 
of three axes. Dislocations lead to a phase mismatch of 
±27r/3 in the frustrated (120°) magnetic order and thus 
attract ±1/3 vortex. Consequently, there are six species 
of dislocations: three directions of Burger's vectors and 
two signs of the vorticity. Dislocations with Burger's vec- 
tors in different directions now interact, the interaction-is 
proportional to the cosine of the angle between them,c3 
cf. Eq. (||^) . The charges are given in table [5[ All species 
of dislocations are equivalent and the final equations are 



dy 2 Jdi = 2(2 -7rJ v )y 2 v , 
dy 2 Jdl = 2(2 - 7T J d - ttJ /9) yl 
dJ v /d£ = ~^y 2 v J 2 v ~ (S/^ypl 
dJ d /dl = -12ir 3 y 2 d Jl 



(52) 
(53) 
(54) 
(55) 



Note that there are no neutral triplets in the triangu- 
lar skyrmion lattice, as opposed to the usual triangular 
lattice. Thus the_first omitted terms are of fourth order 
in the fugacitiesO Eqs. (|5^)-(|55|) differ from the square 
lattice case only in the coefficients. The possible melting 
regimes are thus qualitatively the same. The liquid crys- 
tal phase is hexaticO Bare disclinations do carry vortic- 
ity (±1/3) but this fact is irrelevant for the upper transi- 
tion since the vorticity part of the disclination interaction 
is totally screened by free vortices and dislocations. 

The centered rectangular lattice can be generated from 
the square lattice by tilting the angle 9 between the prim- 
itive lattice vectors away from 9 = n/2 but keeping their 
lengths fixed. This tilting leads to an interaction be- 
tween dislocations with Burger's vectors along different 
primitive vectors. However, all elementary dislocations 
still have the same fugacity and effective interaction with 
their respective antiparticles since they are related by re- 
flection symmetry. The charges are given in table [I]. 



For 9 = 7r/2 we recover the square lattice. Our usual 
ansatz, which assumes that the system can be described 
in terms of screening of q v and qd alone, fails here. We 
cannot reduce the problem to four coupled equations for 
general values of <p but need five: 



dy 2 Jd£ 
dy 2 /d£ 
dJ v /d£ 
dJ±/d£ 



2(2- ir J v )y 2 v , 
2(2-7rJ+/2-7rJ_/2- 
-47r 3 y2 J 2 - 4ir 3 y 2 d Jl 



-8n 3 y 2 d J 2 ± 



with the initial conditions for the J 



MO) 

J±(0) 



(l±coscj>)(3q 2 d (0)/2n. 



(56) 
(57) 
(58) 
(59) 



(60) 
(61) 



Nevertheless there can be at most two transitions (apart 
from disclination unbinding) since there are only two fu- 
gacities. The possible regimes are the same as for the 
square lattice if we replace Jd by ( J+ + J_)/2. Above the 
lower melting transition the system is in a liquid crystal 
phase with quasi-long-range orientational order with re- 
spect to a twofold symmetry. This is a two-dimensional 
nematic phase. Elementary disclinations do not carry 
vorticity so that the upper melting transition is simple. 

The simple rectangular lattice is generated from the 
square lattice by stretching it in the (10) direction. Dis- 
locations with Burger's vectors in the x and y direction, 
respectively, now have different energies, increasing the 
number of independent variables. Dislocations still bind 
±1/2 vortex. There are again N = 5 particles and M = 3 
charges. The final renormalization group equations are 



dy 2 Jd£ 
dy 2 dl /d£ 

dy 2 d2 /dl 
dJ v /d£ 
dJ dl /d£ 
dJ d2 /d£ 



2(2 -nJ v )yl 

2(2-nJ dl -nJ v /4)y 2 dl , 

2(2-TrJ d 2-nJ v /4)y 2 2 , 

-A^y 2 v J 2 v - 2n 3 y 2 dl J 2 - 2n 3 y 2 d2 J 2 ., 

-87T y dl J d i, 

-8tt 3 ^ 2 J| 2 . 



(62) 
(63) 
(64) 
(65) 
(66) 
(67) 



Since the dislocation energies are different, there are two 
distinct dislocation unbinding transitions at T d \ and T d2 
except if they are driven by vortex unbinding. If there 
are two structural transitions, the phase between them 
has translational quasi-long-range order in, say, the x 
direction, but short-range order in the y direction. It 
is thus a two-dimensional smectic phase. There are five 
possible scenarios, which follow from our considerations 
for the square lattice and are not discussed here. For T > 
max(Tdi, T d2 ) we again have a nematic phase (twofold 
rotational symmetry). Elementary disclinations do not 
carry vorticity and the upper melting transition is simple. 

The lower melting transition of the honeycomb lat- 
tice is trivial in this context since its dislocations do not 
carry vorticity: The structural lattice has the same ba- 
sis as the magnetic lattice so that dislocations (with any 
Burger's vector) cannot lead to a phase mismatch. Thus 
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antiferromagnetic quasi-long-range order could persist in 
the liquid crystal phase above the dislocation unbinding, 
transition, which is expected to have unique properties.!^ 
Disclinations carry ±1/2 vortex and any remaining mag- 
netic order is destroyed at the upper melting transition. 

Finally we stress again that these considerations can 
only yield the possible sequence of BKT transitions 
for any lattice type. This approach cannot describe 
other transitions directly, such as first-order melting 
transitionscj or first or second-order structural transi- 
tions. Structural transitions in particular could take 
place since magnetic disordering leads to a reduced ef- 
fective magnetic interaction between skyrmions, affecting 
the stability of the various lattice types in different ways. 

In fact, the question arises of whether the square 
skyrmion lattice can be stable at all above the magnetic 
disordering transition. The answer is affirmative since 
the magnetic part of the interaction is of short range: As 
long as the magnetif-BKT correlation length £„, which 
for T>T r satisfiesEaN 



exp 



(68) 



is much larger than the range of the magnetic interaction, 
£xy ■ the effect of magnetic disordering on the lattice en- 
ergies is negligible. Only when, at a higher temperature, 
becomes comparable to (xy , frustration becomes im- 
portant. In this case we expect the magnetic interaction 
to be effectively reduced so that eventually the triangular 
lattice becomes favorable. If at T = Td still £ v > £xy, 
the square lattice melts and forms a tetratic phase be- 
fore a structural transition takes place. It can even ex- 
hibit the upper melting transition without any structural 
transition taking place, depending on the non-universal 
constant b in Eq. ( p8[ ) . The same kind of argument holds 
for other lattice and liquid crystal structures. 



IV. CONCLUSIONS 

As noted above, the lattice types relevant for present 
quantum Hall systems are the triangular, centered rect- 
angular, and square lattices. A rough estimate of the 
actual core energies and interactions in the real quantum 
Hall system using the model of Sec. II A and results of 
Ref. ~ 



indicates that the low-density, triangular lattice 
typically shows decoupled magnetic and melting transi- 
tions, cf. Fig. ^(a). The magnetic skyrmion interaction, 
and thus the magnetic stiffness and the vortex energies, 
are small since a strong magnetic interaction would make 
the triangular lattice unstable. The effective magnetic 
stiffness is further reduced by frustration. 

On the other hand, the centered rectangular and 
square lattices usually show simultaneous transitions, 
i.e., magnetic order persists up to the lower melting tem- 
perature. Experimentally, the two possible scenarios of 
Figs. ||(b) and (c) are probably not easy to distinguish. 



Note also that vortex driven simultaneous transitions, 
Fig. |J(b) , require fine tuning of magnetic and structural 
stiffnesses so that this scenario is probably rare. 

Quantum fluctuations can affect these results. As 
noted above they are expected to lead to quantum melt- 
ing in the vicinity of the classical square-centered rect- 
angular line. They should also destroy magnetic order in 
the triangular lattice at sufficiently low densities, where 
the magnetic interaction becomes exponentially small. 
More experiments are needed to test the predictions of 
this paper. Sharp structures in the Knight shift, nuclear 
relaxation rates, or specific heat as functions of temper- 
ature would be indications for phase transitions of the 
skyrmion system. Particularly valuable would be exper- 
iments on the electronic susceptibility x(q, uS) as a func- 
tion of temperature and filling factor. 

In conclusion, we have performed a Berezinskii-Koster- 
litz-Thouless renormalization group study of melting and 
magnetic disordering in various lattice geometries in or- 
der to understand the behavior of the skyrmion lattice 
in quantum Hall ferromagnets at finite temperatures. 
The behavior of the skyrmion system is determined by 
the two facts that skyrmions (i) are non-collinear mag- 
netic defects and (ii) carry electrical charge. In the long- 
wavelength limit the in-plane magnetization components 
can be described by a U(l) (XY) degree of freedom as- 
sociated with each skyrmion. The XY "spins" couple 
antifcrromagnetically and can lead to antiferromagnetic 
quasi-long-range order. Dislocations in most skyrmion 
lattice types lead to a mismatch in the XY degree of free- 
dom, which makes the dislocations bind fractional vor- 
tices and leads to coupling of translational and magnetic 
excitations. For most lattice types there are three dis- 
tinct scenarios for the lower melting transition: (i) a BKT 
magnetic disordering transition at a lower temperature 
than BKT melting, (ii) simultaneous transitions where 
the magnetic stiffness shows a universal BKT jump, and 
(iii) simultaneous transitions where the effective disloca- 
tion stiffness shows a universal jump. 

The lattice types we have studied are motivated by 
the possible ground states of a simple classical model 
of the skyrmion system, which uses the large-separation 
limit of their interaction. It shows a surprisingly rich 
T — phase diagram, which suggests that upon increas- 
ing the skyrmion density a frustrated triangular ground 
state first gives way to a centered rectangular lattice with 
Neel order and only at higher density to a square lattice. 
Quantum melting is expected to take place in the vicinity 
of the latter transition. 
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FIG. 2. Classical skyrmion lattice phase diagram at T = 0. 
Thin solid lines are continuous phase transitions, whereas 
heavy lines are first-order transitions. The honeycomb phase 
is denoted by 'h' and the centered rectangular phase by 'cr.' 
The dotted lines in the simple rectangular phase are lines of 
constant anisotropy ij < 1. The employed model is quantita- 
tively correct for [3 2 <C 1. Also, real systems are expected to 
have (3 /a 2 < 1. 
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FIG. 1. Sketch of a dislocation in a square skyrmion lattice. 
The arrows denote the internal U(l) degree of freedom. As 
seen in (a), dislocations lead to a phase mismatch of ±7r in 
the U(l) degree of freedom. In (b) the U(l) angles have been 
allowed to relax and the dislocation has acquired half a vortex. 
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FIG. 3. Dispersion of magnetophonons and XY magnons 
along high-symmetry directions in the magnetic Brillouin 
zone of the square skyrmion lattice. The calculation was done 



for 



0.1 and /3/a 2 = 1. 
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FIG. 4. The three possible scenarios close to the lower 
melting transition of the square skyrmion lattice: (a) decou- 
pled transitions, (b) vortex driven simultaneous transitions, 
and (c) dislocation driven simultaneous transitions, see text. 
The vortex stiffness J v and the effective dislocation stiffness 
Jd + Jv/4: are plotted as functions of temperature. The graphs 
have not been obtained by actual integration but are rather 
sketches meant to emphasize the universal features. 
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TABLE I. Coulomb gas charges of the square skyrmion 
lattice. q v is the vortex charge and q d is the dislocation charge. 



Qv 
q v /3 

-q v /3 
q v /3 

-qv/3 
qv/3 

-qv/3 





qd 
qd 

cos(27r/3)<j<j 
cos(27r/3)g d 
cos(27r/3)g d 
cos(27r/3)<7<i 







sin(27r/3)g<i 
sin(27r/3)g d 

- sin(27r/3)g d 

- sin(27r/3)g<j 



TABLE II. Charges for the triangular lattice. 



qv/2 
-q v /2 

qv/2 
-qv/2 



qd 

cos(6)q d 
cos(8)q d 







sm(Q)q d 
sin(9)q d 



TABLE III. Charges for the centered rectangular lattice. 
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FIG. 5. Sketch of a threefold disclination in the square 
skyrmion lattice. The disclination does not frustrate the an- 
tiferromagnetic order of the XY spins (arrows) . 
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